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Abstract 

QCD sum rules for chiral partners in the open-charm meson sector are presented at 
nonzero baryon net density or temperature. We focus on the differences between pseudo- 
scalar and scalar as well as vector and axial-vector D mesons and derive the corresponding 
Weinberg type sum rules. This allows for the identification of such QCD condensates 
which drive the non-degeneracy of chiral partners in lowest order of the strong coupling 
Us and which therefore may serve as "order parameters" for chiral restoration (or elements 
thereof). 
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I. INTRODUCTION 



Chiral symmetry and its breaking pattern represent important features of strong 
interaction physics. The non-degeneracy of chiral partners of hadrons is considered 
to be a direct hint to the spontaneous chiral symmetry breaking in nature which 
characterizes the QCD vacuum. In fact, the distinct difference of iso- vector- vector 
and iso- vector-axial- vector spectral functions deduced from r decays l|, |2| gives one 
of the empirical and precise evidences for the breaking of chiral symmetry. The low- 
energy strengths of the mentioned spectral functions, concentrated in the resonances 
p(770) and ai(1260), deviate strongly from each other and from perturbative QCD 
predictions. (For a dynamical interpretation of the two spectra see js], Q].) This 
exposes clearly the strong non-perturbative effects governing the low-energy part 
of the hadron spectrum. The spontaneous symmetry breaking is quantified by the 
chiral condensate (gg), which plays an important role in the Gell-Mann-Oakes- 
Renner relation connecting hadronic quantities and quark degrees of freedom (cf. js]. 



Weinberg type chiral sum rules for differences of moments between light vector 



BE 



and for a 



and axial- vector spectral functions have been developed for vacuum 
strongly interacting medium at finite temperature [9^ . It is the aim of the present 
paper to present analog sum rules for the scalar and pseudo-scalar as well as vector 
and axial-vector mesons composed essentially of a heavy and a light quark, e.g. uc 
and dc realized in vacuum as L)±(1867), = 0", D°(1865), = 0", L'*(2400)°, 
JP = 0+, D*(2007)°, = 1-, D*(2010)±, = 1-, Di(2420)°, = 1+ (there 
is no confirmed charged scalar or axial- vector state in the open charm sector) 10 1. 



Such open charm degrees of freedom will be addressed in near future by the CBM 



12j collaborations at FAIR in proton-nucleus and anti-proton- 
Accordingly we are going to analyze the chiral sum rules in 



El- 



ll( and PANDA 
nucleus reactions 
nuclear matter. 

Chiral symmetry is explicitly broken by nonzero quark masses, but quantum 
chromo dynamics is still approximately invariant under transformations which are 
restricted to the light quark sector. The associated conserved currents are the well 
known light quark vector and axial- vector currents. Their mixing under chiral trans- 
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formation and Wigner's non-degeneracy of the hadron spectrum is interpreted as the 
spontaneous breakdown of chiral symmetry. Although currents, related to mesons 
which are represented in the quark model as composed of a light quark and a heavy 
quark, are neither conserved nor associated with a symmetry transformation, vec- 
tor and axial-vector currents still mix under a transformation which is restricted 
to the light quark sector. Consider for example the infinitesimal rotation in flavor 
space i/j e-^*''®"^V ~ (1 - irOT) V', where T e SU(n/), T = 1 for the vector 
transformations and F = 75 for the axial transformations, and denote a set of 
infinitesimal rotation parameters (rotation angles). The vector current transforms 
as 

jj'" = V^Tm^V' ^ ij'" ± ^VS^7M©" r , r]_ , (1.1) 

where -|- (— ) refers to the vector (axial) transformations and r denotes a matrix 
in flavor space. In the 3-flavor sector ip — (■0i, ■02, V's)) where the first two quarks 
are hght, the choice 6 = (61, Os, 83, 0, ... , 0) and 2^ = A", the well-known Gell- 

Mann matrices, clearly leaves the QCD Lagrangian invariant. To be specific, let us 
consider r = (A^ -|- iX^)/2, where the corresponding current transforms as 

jj'" = ^i7m^3 ^ jj'^' = jj'" - 1 {i^ir^.^sQs + (ei + 1Q2) V^2r7;.^3) • (i-2) 

Obviously, the heavy-light vector current mixes with heavy-light axial-vector cur- 
rents if an axial transformation, with F = 75, is applied. Analog expressions hold 
for spin-0 mesons, and the result is the same as in the spin-1 case, but with the 
replacements 7^ — )■ 1 for the vector transformation and i7^[i",T]_ [t"',T]_^ for 
the axial transformation. Hence, a symmetry and its spontaneous breakdown in the 
light quark sector must also be reflected in the spectrum of mesons composed of a 
heavy and a light quark. In particular, the splitting of the spectral densities between 
heavy-light parity partners must be driven by order parameters of spontaneous chiral 
symmetry breaking only. 

Probing the chiral symmetry restoration via the change of order parameters re- 
quires a reliable extraction of their medium dependence. As is long known for the 
vacuum case, the chiral condensate is numerically suppressed in QCD sum rules in 
the light quark sector due to the tiny light quark mass, but occurs amplified by the 
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large heavy quark mass in QCD sum rules involving a light and a heavy quark [l^ . 
Despite of the amplification in terms of mduu), with rric and (uu) to be evolved to 
the appropriate scale, in case of the D mesons, the dependence of the in-medium 
D meson spectrum on the chiral condensate is not as direct as anticipated. This is 
clear in so far as there are always particle and antiparticle contributions to the spec- 
trum of pseudo-scalar states and one has to deal with both pieces. This accounts 
for inherent suppressions and amplifications of different types of condensates due to 
the generic structure of the sum rule in that case. Indeed, a precise analytic and 
numerical investigation points to competitive numerical impacts of various conden- 
sates for the mass splitting of particle and antiparticle. Whereas, the determination 
of the mass center rather depends on the modeling of the continuum threshold 15 1. 

In the difference of chiral partner spectra (Weinberg type sum rules) the de- 
pendence on chirally symmetric condensates drops out. At the same time, the 
amplification of the chiral condensate by the heavy quark mass is still present. This 
makes Weinberg type sum rules of mesons composed of heavy and light quarks an 
interesting object of investigation. 

Our paper is organized as follows. The operator product expansion (OPE) and 
a suitable projection are elaborated in section [III for differences of current-current 
correlators. These are spelled out for chiral partners of open charm mesons in the 
form of Weinberg type sum rules with the focus on the OPE sides in section llllj 
a remark on heavy-quark symmetry is included as well. Numerical examples are 
exhibited in section IIII D[ The summary can be found in section IIVI In Appendix 
|A| some relations between current-current correlators are considered. Appendices 
IB] and O list projection coefficients and prove a necessary relation for the quark 
propagator. 



4 



II. DIFFERENCES OF CURRENT-CURRENT CORRELATORS AND 
THEIR OPERATOR PRODUCT EXPANSION 



We consider the currents 

f{x):=qi{x)q2{x), (2.1a) 

fix) := qi{x) q2ix) , (2.1b) 

j^{x) ■■= qi{x) 7m q2{x) , (2.1c) 

j^ix) ■■= qi{x) 757m q2{x) (2. Id) 

and the corresponding causal correlators 

n(S'P)(g) = z j d^xe^^^iT [f^'''Hx)f^'^'>^{0)]) , (2.2a) 

^Z'^\q) = ^ j d^a;e^«^(T [jf '^^^0)] ) , (2.2b) 

where T[. . .] denotes time-ordering and (. . .) means Gibbs averaging. Retarded and 
advanced correlators are related via analytic continuation to the causal correlator 

Q. 

In the rest frame of the nuclear medium, i.e. n = (1,0), and for mesons at rest, 
i.e. q = (1,0), the (axial-) vector correlator can be decomposed as 

n^'^n^) = - 9,.) 4^'^)(g) + ^nr)(g) (2.3) 

with U^^'^\q) = I - (7^") niy)(g) and Il^T^'^\q) = ^q''q'Il'll'^\q) written 

explicitly covariant. Ut contains the information about (axial-) vector degrees of 
freedom. Ul refers to (pseudo-) scalar states and can be related to the (pseudo-) 
scalar correlator (see Appendix A for details). 

We now proceed with the OPE's for G jn^^.P), n^^'^) = c/'^^nS^'^^ } . We do 
not include the as corrections which would arise from inserting the next-to-leading 
order interaction term in the time-ordered product of the current-current correlator. 
Such terms account, e.g., for four-quark condensates. They are of mass dimension 6 
and beyond the scope of this investigation. According to standard OPE techniques 
19|). the time-ordered product can be expanded into normal-ordered 



(see e.g. 
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products multiplied by Wick-contracted quark-field operators. Dirac indices can 





(a) (b) ^ ^ (c) /^~^ ' 

2 2 2 

FIG. 1: Diagrammatical representation of (a) n^'^^ and (b) n^^^. (c) is for U^^^ with two 
gluon lines attached to the free propagator of quark 1. Double lines stand for the complete 
perturbative series of the quark propagator, wiggly lines are for gluons, and circles denote 
non-local quark condensates, while crosses symbolize local quark or gluon condensates. 

be projected onto elements of the Clifford algebra, T G {1, 7/^, o"^;/, 757/^, 75}, which 
provides an orthonormal basis in the space of 4 x 4 matrices with the scalar product 
{A, B) = ^Ttj:,[AB]. Color indices can be projected onto an analogously appropriate 
basis. Thereby, color and Dirac traces of the quark propagator occur. Using the 



background field method in fixed-point gauge (cf. [l7|, [iSj]) for the gluon background 
field Afj^{x), x^A'^^x) = 0, one can expand the non-local quark field operators by a 
covariant expansion q{x) = ^"^nf (-^ai • • • ^ar,Q)\x=o arrive, after a Fourier 
transformation, at 

1 °° 1 
r n=0 

+ (-l)"g2rTrD [rr^5i(g)r^] D^^q^ :) 
= U^^'^\q) + U^^^\q) , (2.4) 

where Z)^,^ = . . . Da„ (with an analog notation for the partial derivative) and 
quark fields and their derivatives are taken at x = 0. For X denoting vector 
and axial- vector states, (r^^).^. (r^^)^^ = (F;^).,^. (r^^)^, is understood. n(°)(g) de- 
notes the fully contracted (depicted in Fig. [T] (a)) and n(^)(g) the 2-quark term 
(see Fig. [T] (b) for {tjj . . . ip)) , Trc,D means trace w.r.t. color and Dirac indices. 
In doing so, we introduce the momentum-space perturbative quark propagator in 
a gluonic background field as S{q) = f d^x e^'^^S{x,0) = J d^x e~^'^^S{0,x) and 
S{x,y) = — i(0|T [q{x)q{y)] |0). To have the same structures in both quantities U^^^ 



and n(^), we also project the matrix product F^ST^ onto this basis (see Appendix 

n(P.s)(0)(^) = _,J ^(i^Trc [TT^[S,ip + q)]TT^[S,ip)] (2.5a) 

+ (- 1) (P'S)TrD [^2 {p + qh,]TTB [Si {p)Y] 
+ ^Ttb[S2{p + g)a^.]TrD[Si(p)an 
+ (- 1) (P'S)TrD [^2 {p + g)757M] Trc [Si [pHY] 
+TrD[52(p + g)75]TrD[5i(p)75]]:), 
n(v,A)(o)(^) = , I ^(:Trc [Tri,[S2ip + q)]TrB[Siip)] (2.5b) 

+ (-l)(^'^)TrD[S2(p + g)7M]TrD[Si(p)7^1 

- (-l)(^'^)TrD[52(p + g)757M]TrD[^i(p)757l 
-TrD[52(p + g)75]TrD[5i(p)75]]:) , 

where (—1)*^^'^-' = —1 for pseudo-scalar and vector mesons and (— 1)(^'^) = 1 for 
scalar and axial- vector mesons. 11"^*^^) (g) in (12. 4p may be simplified using Tab. HTlin 
Appendix [Bl 

To obtain the OPE's for the difference of chiral partners n^~^ = 11^ — n"^ 
and n^~^ = — 11"^, one can use (12. 5p or directly project the occurring anti- 
commutators. The result reads 

nP-sW(g) = -z I ^(:iTrc{TrD[52(p + g)]TrD[5i(p)] 

+ ^TrD[52(j9 + q)a,,]TTr>[Si{p)a^n +'Ttb[S2{p + qH]TTi,[Si{ph,]}:) , 

(2.6a) 

np-s(2)(,) = ^ i_!L^ ^ (%^.^„ra- (TrD[r52(g)]) qi 

n r 

+ q.Td''- (TrD[r5i(-g)]) '^^^q^:) (2.6b) 

n^-^(°n?) = ^ / (^(:2Trc{TrD[52(p + g)]TrD[5i(p)] 

-TrD[5'2(p + g)75]TrD[5i(p)75]}:) , (2.6c) 
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n^-A(2)(^) = _^Li2 J2 x(:giK^„r9"^"(TrD[r52(g)])gi 

n r 

+ q,rd^- (TrD[r5i(-g)]) D,^„g2:) , (2.6d) 

where the Chfford basis is modified now by the imaginary unit in front of 75 for the 
vector-axial-vector difference. The advantage is that we are left with three different 
types of Dirac traces for the quark propagators in the P— S case and only two in the 
V— A case. 

The perturbative series for a momentum-space quark-propagator in a gluon back- 
ground field in fixed-point gauge is presented in Appendix O Its derivative is given 
by the Ward identity d^S{q) = —S{q)T^{q,q;0)S{q), where T^{q,q;0) denotes the 
exact quark-gluon vertex function at vanishing momentum transfer. r^(g, q; 0) = 7^ 
holds for a classical background field meaning that the Ward identity for the com- 



20| 



plete perturbation series has the same form as for free quarks 

For the limit of a massless quark fiavor attributed to qi, mi — )■ 0, one can show 
(see Appendix O that TrD[r 5*1(5)] = for F G {l,cr^;^,75} and 7^ 0. In this 
limit only the diagram in Fig. [1] (b) gives a contribution to 11'^^^ The corresponding 
diagram for 1 < — > 2 vanishes, because the sum over Dirac matrices in (12. 6p covers 
such elements where the corresponding traces vanish. Hence, for chiral partner sum 
rules the often used approximation of a static quark, which results in vanishing 
heavy-quark condensates, is not necessary as their Wilson coefficients vanish in the 
limit of the other quark being massless. This means that, if one is seeking the 
occurrence of quark condensates in lowest (zeroth) order of the strong coupling 
as, at least one quark must have a nonzero mass. (Higher order interaction term 



insertions cause t 
of a,, cf. 



m 



le occurrence of further quark condensates proportional to powers 



22j for examples.) Hence, the structure of the OPE side of the 
famous Weinberg sum rules for light quarks w.r.t. quark condensates is shown in all 
orders of the quark propagators and quark fields. 

For the completely contracted term H*^*^^ the situation is somewhat more involved. 
In case of two light quarks, i.e. mi 2 — ?■ 0, the limiting procedure and the momentum 
integration commute and, hence, H*^*^^ = is obvious, because all non- vanishing terms 
drop out in the chiral difference. If one of the quarks has a nonzero mass, the limiting 
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procedure and the momentum integration do not commute due to the occurrence 



of infrared divergences 



23I-27 



-27|, e.g. for the term depicted in Fig. [1] (c) which is 
proportional to the gluon condensate. As the integration domains in fl2.5p and fl2.6a[ 
I2.6cp involve momenta p = 0, Tr[r5'(p, m)] does not converge uniformly for m — )■ 0. 
Hence, the integration and the limit m — )■ cannot be interchanged. In fact, a direct 
calculation of the perturbative contribution to chiral partner OPE's shows that the 
infrared divergences are the only remaining terms. All finite terms cancel out in the 
chiral difference. These divergences have to be absorbed by introducing condensates 
which are not normal ordered. We demonstrate here explicitly their cancellation. 
Up to mass dimension 5, the only product of traces which contributes to 11*^°^ is 
TrD[5'i]TrD[5'2]. This is true for arbitrary quark masses. Indeed, up to order a] 
one can show that the contribution — J d'^pTTY)[S^\p + q)]{:TTB[S^\p)]:) , which is 
the remaining term of (12. 5p after the limit mi — )■ has been taken, is canceled by 
/ d'^pTT-D[S^\q)]{TT-D[S^\p)]) . The latter quantity is introduced by the definition 



of non-normal ordered condensates 

{Q\qd[D^]q\Q)={Q\ ■.qOlD.jq: \Q) 

d^ 



J ^(n|Trc,D [O [-tp, - tA,] S,{p)] \n) 



(2.7) 



where O [D^] is an arbitrary operator depending on D^, and denotes the Fourier 
transform of A^, resulting in a factor J {Tty)[S^\p)]) to the Wilson coefficient 
TrD[5'2°^(g)] of the chiral condensate in (12. 4p . Hence, in case of heavy-light systems, 
first the integration has to be performed, then one has to introduce non-normal or- 
dered condensates according to (12. 7p and afterward the limit mi — )■ can be taken. 
In case of equal masses, mi = m2, the divergences cancel each other with virtue to 

J d'p [Trn[S['\p + q)]Trn[s!t\p)] + Ttb[S['\p + q)]TTj,[Si'\p)]) = . (2.8) 

If two heavy (static) quarks are considered, only H*^") gives a contribution to the 
chiral OPE, whereas for two massless quarks, both H^^^ and H^^^ vanish. 

Putting everything together this means that for light quarks (mi^2 ^ ^qcd) in 
the chiral-difference OPE the corresponding traces and, therefore, the corresponding 
Wilson coefficients vanish, while for heavy quarks (mi_2 ^ ^qcd) the condensates 
vanish. To obtain quark condensates in order the two flavors must be of different 



mass scales, i.e. qi G {u,d{,s)} is a light quark and q2 G {c,b,t} is a heavy quark. 
Hence, to seek for condensates which are connected to chiral symmetry breaking as 
possible order parameters in order a natural choice is to consider chiral partner 
mesons composed of a light and a heavy quark. Given the above mentioned ex- 



perimentally envisaged research programs at FAIR |llHl3| we focus on open charm 
mesons. The presented formulas may be directly transferred to open bottom mesons 
by rric mi,. 

III. CHIRAL PARTNERS OF OPEN CHARM MESONS 

We now consider a light (gi = q) and a heavy (g2 = qc) quark entering the 
currents in Eq. (12. ip . 

A. The case of P-S 

For the P— S case we consider the pseudo-scalar D{0~) (-D^, and Z)°) and 
its scalar partners Dq{0~^). Of course, all the results also account for other heavy- 
light (pseudo-) scalar mesons. (For open charm mesons Dg which contain a strange 
quark, however, the limit — )■ may not be a good approximation and terms 
oc rris should be taken into account as well.) 

Up to and including mass dimension 5, after absorbing the divergences in non- 
normal ordered condensates, the OPE gets the following compact form 

nP-S(g) = Ifiq) - nS(g) = U^^^\q) - U''(^^\q) 

. (3.1) 

n ' r 

where the sum over the elements of the Clifford algebra does not contain 75 up to 
this mass dimension anymore. To evaluate the condensates in (13. ip in terms of ex- 
pectation values of scalar operators, Lorentz indices have to be projected onto g^y 
and e^ysA in vacuum, whereas the medium four-velocity provides an additional 



structure at finite densities and/or temperatures [28|]. Hence, new condensates must 
be introduced which vanish in the vacuum. Thereby, temperature and density de- 
pendences stem from Gibbs averages of medium specific operators. The evaluation 
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in the nuclear matter rest frame = (1,0) for mesons at rest yields 



qi - ml {ql - mlY 

+ [{qga^q) - 8{qDlq)] J^^^ (3-2) 

where we separated a medium-specific term (last hne, (qga^q) — 8{qDlq) = (A)) 
vanishing in vacuum, g'^ = Anas, cf^^ = i [7^, ■ju] /2, = i [-D^, D,j] and q denotes 
either d or u quark field operators. The Eq. (13.21) reduces to the vacuum result 
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29 



30[ at zero density and temperature. The condensates (qq), {qga'^q) and 
(qD^q) may have different medium dependences. 

An odd part of the OPE, i.e. U{qo) = ff™^(go) + qoU"'^'^ (qo) with n^^^"(go) = 
(n(go) + n(-go))/2 and W'^'^iqo) = {U{qo) - n(-go))/(2go), where both functions 
are even w.r.t. go, does not appear up to this mass dimension. Although there is no 
7^ projection of the condensates for difference OPE's of chiral partner, it may arise 
from an odd number of derivatives in f l3.1l) . 



As usual within a QCD sum rule analysis |17H19|. a Wick rotation to Euclidean 
momenta, go = has to be performed. The OPE in the deep space-like region 
based on asymptotically free quark degrees of freedom can be expected to con- 
verge. The OPE of the current- current correlation function in turn is related to 
the hadronic spectral density by an in-medium dispersion relation. To evaluate 
the correlator at arbitrary values of the energy go off the real axis, Imgo 7^ 0, e.g. 
go = —Q^ —7- —00, an integral over its discontinuities (spectral densities) along 
the real axis An(go, ^ = 0) = ^ lim^^Q [n(go + ie) — n(go — i^)] = Inin(go) is per- 
formed. After a Borel transformation the result is 



-2mc(gg) 



(3.3) 



M is the Borel mass. In a medium, the integral over the spectral densities runs over 
positive and negative frequencies u, covering both particles and their antiparticles, 
of given quantum numbers. 

It is instructive to cast Eq. (13. 3 P in the form of Weinberg type sum rules 
This can be accomplished by expanding the exponential on both sides and comparing 
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the coefficients of inverse powers of the Borel mass. In such a way we can relate 
moments of the spectral P— S differences to condensates via 



1 

vr 
1 
7r 
1 
vr 



du ujAU^-^{uj) = - 2m^{qq) , (3.4a) 

— oo 
oo 

duLU^Alf~^{u) = - 2ml{qq) + m^{qga^q) - (A) (3.4b) 



— oo 
oo 



duJUj''AU^-^{uj) = - 2ml{qq) + 3m'^{qgaWq) - Sm"^ {A) + . . . , (3.4c) 



where the dots denote the neglected contribution of mass dimension 7. This gener- 
alizes the OPE side of Weinberg type sum rules to scalar and pseudo-scalar mesons 
in the heavy-light quark sector for the first time. 

If one attributes chiral symmetry to the degeneracy of chiral partners (i.e. the 
l.h.s. of (13. 4 p vanishes) the vanishing of {qq) and {qD^q) = | {{qga'^^q) — (A)) on 
the r.h.s. is required. In this spirit, these condensates may be considered as possible 
order parameters of chiral symmetry. Note that (13.41) also allows to consider the 
omitted mass dimension 7 condensate as an order parameter. Remarkably, the 
chiral condensate (qq) of light quarks figures here in conjunction with the heavy 
quark mass as parameter for chiral symmetry breaking in each of the moments (for 
vacuum, cf. [29]). Of next importance is (qD^q), again in combination with the 
heavy quark mass. The r.h.s. quantities must be taken at a proper renormalization 
scale. Formally, in the chiral (i.e. strictly massless) limit for all quarks the r.h.s. of 
(13. 4 p would vanish. 



B. The case V-A 



In the same manner we proceed in the V— A case. From Eqs. (I2.6cp and (]2.6dp 
we obtain up to and including mass dimension 5 

B^-^q) = nV(g) - n^(g) = U''^^\q) - U'-^'^q) 

{1^5} ^ ^ (3.5) 

n' 



(_-\n {I'^-^s} 

EM^2 J2 (giK^„r9-TrD[r52(g)]gi) , 



n 



where only the 1-projection survives. The in-medium evaluation results in 
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and together with f lA14p and (13.20 we obtain for the correlator containing the infor- 
mation about the vector and axial-vector degrees of freedom 



V-A 



(g) = n^-s + 



rrir 



(g2 -m2)2 



{qgaGq) + - 



rrir 



3(g, 



m 



2\2 



(A). 



(3.7) 



For vacuum, the result of 
given by 



14 



30| is recovered. The Borel transformed sum rule is 



-2mc{qq) 



171" 



2M4 



(qga'^^q) - (A) 



4 rrir 



2M4 3 M2 



(3.8) 



The corresponding moments therefore are 

(w) = - 2mc{qq) 



1 

TT 
1 
71 
1 
TT 



oo 

hOO 



2ml{qq) mc{A) 



— oo 
oo 



duju^AU}(.-^{uj) 



— 2ml{qq) + m 
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(3.9a) 
(3.9b) 

ml{A) + ... (3.9c) 



with the same meaning of " . . ." as above. This second set of Weinberg type sum 
rules contains the same condensates as the first set in Eq. (13.41) but in different 
combinations. Again, addressing chiral symmetry to the l.h.s. one may consider 
the chiral condensate (gg) and (A) as possible order parameters. As the P— S case 
allowed us to identify {qD^q) as order parameter, {qgaWq) also qualifies as order 
parameter. 
In addition. 



dcJwAIl, 



V-A/ 



[UJ 







(3.10) 



follows from ( 12. 3p but for the decomposition II^j^ = {q^qi, — q'^gfj,u)IlT + qfiqu^h- This 
corresponds to Weinberg's first sum rule f?]. Note that, in contrast to Weinberg's 
original sum rule I?], no Goldstone boson properties appear on the right hand side 
of (I3.10p because the heavy-light currents involved in our case are generally not 
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conserved. The Borel transformed sum rule for IIt = Ht / reads 



rrir. 



1 _ (^-ml/M^ 
1 



{qga'^q) - 



1 _ e--^l/M- I 1 



+ (A)- 



c . 

^2 /A,f2 



M2 2MV_ 



+ 



2M4 



(3.11) 



and reproduces ( 13. 9p . Note that using IIt, instead of IIt, is more appropriate for 
the heavy-quark hmit. 



C. Heavy-quark symmetry 



In the heavy-quark hmit, m| ^ the leading contributions in (13. 7p are 



m|>|i}2| 



{qq) 

1712 



(3.12) 



where o 'i 



/ 2 \ 

~^S^o (^) been exploited. This result is in agreement 
with the expected degeneracy of vector and pseudo-scalar mesons and of axial- vector 
and scalar mesons in the heavy-quark limit 3l|- Actually the (gg) parts of 11^^^ 
and n!^~^ agree as one can check by comparing (13. 3p with (13. 8p . 



D. Numerical examples 



Using the condensates from [15] the r.h.s of Eqs. ( 13. 2p . (13. 8p and (13. lip , i.e. the 
OPE sides, are exhibited in Fig. [2] as a function of the Borel mass for vacuum and 
cold nuclear matter. The Borel curves are significantly modified by changes of the 
entering condensates due to their density dependences, i.e. their magnitudes are 
lowered by approximately 30% of their vacuum values. It should be emphasized, 
however, that the density dependence is estimated in linear approximation. For the 
chiral condensate {qq) it is known 32|, |33| that the actual density dependence at 
zero temperature is not so strong. One may expect accordingly a somewhat weaker 
impact of the medium effects. 

For further evaluations, the hadronic spectral functions or moments thereof must 



be specified, e.g. by using suitable moments as in 



34 



35| 
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M [ GeV ] 




M [ GeV ] 




FIG. 2: The OPE sides of Eqs. 1^ (upper panel, P-S), (middle panel, V-A) and 
(j3.1ip (lower panel V— A) as a function of the Borel mass in vacuum (solid red line) and at 
nuclear saturation density n = 0.15 fm^'^ (dashed blue line). The employed condensates 
are {qq) = (gg)o + 45n/ll, {qga'^q) = 0.8 GeY^ {qq) o + 3n GeV^, {qq)o = (-.245 GeV)^ 
and (A) = 8 X 0.3 n GeV^. 
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IV. SUMMARY 



In summary, we present difference QCD sum rules for chiral partners of mesons 
with a simple quark structure. Focussing on the operator product expansion of the 
current current correlator in lowest (zeroth) order of an interaction term insertion, 
only the combination of a light and a heavy (massive) quark yields a non-trivial 
result: Differences of spectral moments between pseudo-scalar and scalar as well 
as vector and axial-vector mesons for condensates up to mass dimension < 5 are 
determined by the combinations mc(gg), mc{qg(J'^q) , and mc{{qg<7'^q) — 8{qDlq)) 
(to be taken at an appropriate scale) which may be considered as elements of "order 
parameters" of chiral symmetry breaking (see 
mixed quark-gluon condensate at finite temperature). Vanishing of these conden- 
sates at high baryon density and/or temperature would mean chiral restoration, i.e. 
the degeneracy of spectral moments of the considered chiral partners. Chiral part- 
ners of mesons with light-light or heavy-heavy quark currents are non-degenerate in 
higher orders of a^, as exemplified by the Kapusta-Shuryak sum rule. The famous 



361 ] for a lattice evaluation of the 



Weinberg sum rules generalized to a hot medium in Iq] refer to lower mass dimension 
moments where the OPE side vanishes in the chiral limit. 

Our results show a significant change of the OPE side when changing from vacuum 
to normal nuclear matter density. This implies that also the hadronic spectral 
functions may experience a significant reshaping in the medium. 

The open charm meson sector will be investigated in precision experiments by 
two collaborations in near future at FAIR in proton and anti-proton reactions at 
nuclei as well as in heavy-ion collisions. These experiments will shed some light on 
the mass spectrum of open charm mesons and its medium modifications. 

Acknowledgments: The work is supported by GSI-FE and BMBF 06DR9059. 

Appendix A: Relations between current- current correlators 

We denote the currents in Eq. (12. ip by 

jX,r ^^Y'^T^, G {1,^75,7/., 757/.} • (Al) 
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where = collects the fields of n/ quark fiavors and r is an arbi- 

trary matrix in flavor space. Due to the special choice of the current fulfills 
= j^''^\x) = ipT^T'^ip. Lorentz indices are suppressed here for V and A 

currents. 

The currents and are supposed to carry the quantum numbers of scalar 
and pseudo-scalar mesons, respectively. The issue is more complicated for and 
j^. To see this we consider the causal correlator of two vector-field valued, mesonic 
(bosonic) Heisenberg operators and B^, (to be identified with j^''^^''^{x) and 

n^.(g) = ^ j d'xe'^^{T [A^ix)B,iO)]) . (A2) 

A system at nonzero temperature and/or (baryon) density is characterized by a 
heat-bath vector (encoded in the Gibb's averages in (12. 2p ) which can be normal- 
ized by = 1. The heat bath is at rest in a reference frame with n = 0. For a 
system with a conserved current j^, e.g., the baryon current, one can use oc jj^ 
(Eckart choice). For a thermal system without conserved currents one must lock the 
fiow with local energy density (Landau choice, cf. [37]). 

A rank-2 tensor in four dimensional Minkowski space-time can be decomposed 
into six independent algebraic invariants A, . . . ,F being functions of g^, n ■ q: 



n^.(g) = {g,. - ^) a + q^q^B + nlnlC 



+ {q,nl + nlq,) D + {q^nl - n^q,) E + e^J^^q^n^ F (A3) 

where '■= — ^ q^. From invariance w.r.t. parity one obtains F = 0. Imposing 
aAdditionally, time reversal invariance together with translational invariance of the 
considered medium, e.g. nuclear matter, give rise to the symmetry property II^j, = 
Il.u^l1 hence, E = F = {). If the correlator (causal, advanced, retarded correlator 
or spectral density) was four-transversal in both indices, i.e. g'^II^j, = q^H^u = 0, 
E = would follow without additional symmetry constraints. Current conservation 
alone is not a sufficient requirement for the transversality of the causal correlator. 

The coefficients A and C carry the information about the four-transversal degrees 
of freedom referring to vector or axial- vector components which can be decomposed 
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into three-momentum transversal and three-momentum longitudinal states [38|] . The 
coefficient B is related to the four-longitudinal states referring to a scalar or pseudo- 
scalar component. The coefficient D encodes the mixing between three-momentum 
longitudinal (axial-) vector states and (pseudo-) scalar states occurring due to broken 
rotational invariance for an excitation (hadron) moving with nonzero velocity in the 



medium 
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40|. Thus, the spin is no longer a conserved quantum number. 



In the special case n = (1,0) and q = {qo,0), one has = and, therefore 



C and D drop out. Hence, the correlator can 



De decomposed into a longitudinal 



(IIt = —A) and a transverse (IIl = q'^B) part IJ], see Eq. fl2.3p . At zero density and 
zero temperature IIl would vanish for B^, being conserved as demonstrated below. 

To relate the four-longitudinal part of the (axial-) vector current to the (pseudo-) 
scalar current let us consider the correlator 

^>.uiq,p) = j d^xd^ye^«^e-^^'^(T [A^{x)BM]) ■ (A4) 
Translation invariance allows to factor out a momentum dependence 

n^.(g,p) = ^(27r)^5(^)(p - q)Ii,M ■ (A5) 

Contracting with q^^Pu-, performing a suitable number of partial integrations in each 
term and using d^Ao{x) = 9^A^(x)— 9^ylj(x) gives the non-anomalous Ward identity 
relating the longitudinal part of the causal correlator of two Lorentz vectors to the 
correlator of their divergences 



g>^n^,(g,p) = -t' j d^xd^i/e*^^e-^^^(T [zd^A,{xWyB,{y)] 

+ ^S{xo-yo) (q^[A,{x),Bo{y)l^^^^+[Ao{x),zd^^B^{y)]]). (A6) 



The ffist and third terms on the right hand side of this equation vanish if d'^B,y = 
holds. 

Identifying and B^, with the currents above, j^''^(x) := ip^x) '~^^t ip{x), 
j^'^{x) := ijj{x) 757^r^/'(x), and using the equations of motion for the quark fields, 
their divergences are given by 

,5.^-(V,A),. = ^ (^^) [^^ M]^ ^ = (_z)(S,P)/S,P),[.,M]^ ^ ^A7) 
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where the upper (lower) sign corresponds to the vector (axial-vector) current and 
M = diag (mi, . . . , m„^) is the quark mass matrix. 

The equal-time commutators in (1A6P can be obtained from the field anti- 
commutators or, equivalently, from current algebra considerations for the first com- 
mutator as the second involves a divergence. The equal-time current commutation 
relations are given by 



^JV,A),r' 



(y) 



V,[r,r'1-] 



X] 



(A8a) 
(A8b) 



We here consider commutators of currents with the same parity. Setting A^{x) = 
j^^'^'''^(x) and By{y) = j!)^'^^''^\y) in flA6p . inserting the equal time commuta- 
tors (lASp . using translational invariance and comparing with (lASP we obtain the 
non-anomalous Ward identity which relates the causal correlator of (axial-) vector 
currents to the causal correlator of their divergences 



n(s-p) 



r, r 



q 



r.M|,.|r'.M|')-(,-,J-l-"W(/'f"'""'lJ-). 



(A9) 



We turn now to an n/ = 2 fiavor system. Let heavy-light meson currents be 
given by r = r' = (a^ + ia'^)/2, where a* = 2t* are the Pauli matrices, which 
gives fl2.ip . The diagonal mass matrix can be written as M = {Cycr^ + Ca1)/2 



r, 



and the commutators fulfill [r, M]^ = tC{y^a)T, [t, t^-^ = t^^'^\ 
=FC*(v,a)T*-^'^-', where we have defined Cy = mi — m2, Ca = mi + m2 and = a^, 
= 1. Note that qi is attributed to a light-quark field (e.g. up or down) and q2 
to a heavy-quark field (e.g. charm or bottom). 

With these relations the longitudinal part of (12.31) is given as 

g2n[^'^)(g) = q'q'^^Z'^Kq) = Cfv,A)n(^'^)(g)-(#a3^)+qv,A)(^r(V'A)^) (AlO) 



with q = 7^g^. (Equation flAlOp also holds in a three-fiavor system with two light 
quarks and one massive quark.) Furthermore IIl is zero if the current j^'''^' (y) (or B,^ 
in (]A2p ) is conserved, i.e. mi = m2 for vector currents or mi+m2 = for axial-vector 
currents (first and third term on the right hand side of Eq. (lAlOp vanish), and the 
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difference of light and heavy net quark currents being zero, i.e. {'ipqa^ip) — 0. Current 
conservation alone is not sufficient for a vanishing longitudinal projection. Indeed, 
the transversality in case of current conservation generally reads g'^H^Y'^^ {q |t, r') = 
~{ju (^0 ''"'^] _ ) )' which only vanishes in a medium which is symmetric w.r.t. 
the quark flavors of the meson current. If IIl is zero, the transversal projection 
U.^'^\q) is proportional to the trace of the correlator 11^'"^^ (g) = —g'^''Il^f^''^\q)/3. 
Otherwise the trace 11^^'^^ = g'^''U.^^'^\q) contains pieces of (axial-) vectors and 
(pseudo-) scalars. The pure (axial-) vector information is encoded in U.!^'^\q) for 
which one obtains 

3n5^,A) ^ ^^iS,F) ^ ^(^^^3^^ ^ ^^^^(V,A)^^ _ g,u^(y,A) (AH) 

which relates 11^'^^ to the trace of the correlator and 11^^'^^. For mi and m2 being 
arbitrary quark masses, the chiral difference n^~^(g) and the sum nS+P(g) enter 
the chiral difference q^q^ll^i^^iq) and, hence, n!^~^. The expectation value of the 
quark current cancels out in any case: 

nV-A = _1 + ^i j^P-S ^ o ^l^2 j^p+s 

^ 3 \ g^ q'^ 

2 n^-^\ 

+^ {mi{q2q2) + m2{qiqi)) + — ^ . (A12) 
If one quark mass is zero, mi — > 0, one obtains the relation 

q^q^'Ul'^q) = -mjU^'-^iq) - 2m2{mi) (A13) 
and therefore 



n^-^g) = -0n-«(,) - \n^-\q) - ^^{M.) ■ (Ai4) 



The last term contains the interesting combination of light-quark condensate and 
heavy- quark mass. 

Appendix B: Dirac projection 

Consider the trace over Dirac indices in TrD[-S'275{-S'i, 75}]. To write this as a 
product of two simpler traces, one can project the commutator {S, 75} = 
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^6 

yj 


{^,75} 


{7^57^,75} 








iTrn [575 1 


-2TrDf51 

^ L J 


















Ai^<^ 


iTrD[575a^<1 









A\ 












^5 




2TVd[5] 






TABLE I: Coefficients of the projection of {5,75}. 






1 


75 


7/. 


757aj 








TrD[5] 


-TrD[5] 




iTrD[57,,] 


-iTrD[S'7!.] 


-iTrD[57,,] 


-iTrn \S-i,A 
















iTrD[5'757;,] 


-iTVD[5757A.] 


iTrD[5757^] 


iTrD[S'757A.] 




iTrD[S'75] 


iTrD[575] 


-TrD[575] 


-TVd[575] 



TABLE IL Coefficients of the projection of V\^SV\ 



with coefficients Aq = ^Tt-d[{S, 75)0], O being an element of the Chfford algebra. A 
similar expression applies for TrD[5'275{7^S'i7^, 75}] with coefficients listed in Tab. [H 
Alternatively, the coefficients for projecting F^ST^ = Y2d given in Tab. [TTl 

Eq. (12. 5 p can then be obtained by using 750"'^'' = je^uapO''^'^ , leading to 
TrD[5i75a^.]TrD[5275^^1 = TTj,[S,a^,,]TTj,[S2an 

Appendix C: The perturbative quark propagator 

The perturbative quark propagator in momentum space in a weak (classical) 
gluonic background field can be written as S{p) = J2n=oS^''Kp) with S^'^^p) = 
(-l)"5("-i)(p) (^7i) 5(0) (p) = (-l)'^S(o)(p) (^7i) 5("-i)(p) . A denotes a derivative 
operator which arises due to the Fourier transform of the perturbation series for 
the quark propagator in coordinate space from the gluonic background field A^; 
4 = E^=o 4"^ with = {D^, . . . D^„%M) d'^d-^ . . . d--. The trace 
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of each term S^'^\p) can be written as 

TrD[r5W(p)] = (-1)" (\) • • • (^"^.^^- 



X 5^ TrD[r„r^(°np)7^l((5'^"5"^^"TrD[r„-ir„5(o)(j9)7'^"-^]... (CI) 

(9"i9"^iTrD[ri5(°Hp)])^ 



ri,...,r„ 

X 

The sum runs over elements of the Chfford basis. For the sake of a concise notation 



x=0 



and 9°'^'= = . . .9"^ 



we have defined -Da-;%i. = -9 \.i(^l+2) ( -^"i • • • -^"/^mj' 
With these reduction formulas it is only necessary to consider traces of the free 
propagator S^^^ = (p + Tn)l{jP' — w?). Moreover, the single traces form a chain 
making the sum over Fj dependent on traces surviving from the sum over Fj_i. 
In the limit m — )■ and 7^ we have TrD[FiS'*-°''(p, m = 0)] = VF 7^ 7^. 
Hence, from the last trace in (]C1|) only 7^ is passed to the next trace. For the 
latter one we therefore have TrD[FiF2S'(°^(p, m = 0)7^] = VF2 ^ {7a,757o}- As 
TrD[757aF„S'(°)(p,m = 0)7^] = VF„ ^ {7a>757a} also holds, each sum merely 
covers 7^ and 757^. Finally, we obtain for the first trace in the second line of f lCip 
Tr[F„FS'(°^(p, m = 0)7^] = VF ^ {7o,757a}, which means that TrD[F5''^"^(p, m = 
0)] = VF ^ {7M)757At}- As this is true for all orders of the perturbative sum, we 
can conclude that 1i-o\rS{p,m = 0)] = VF ^ {7^,757^}- 
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